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Saudi  Arabia  is  enriched  by  many  geothermal  resources  that  are  located  along  the  western  coastal  part  of 
the  Red  Sea,  in  the  form  of  a  number  of  hot  springs  and  many  volcanic  eruptions.  Wadi  Al-Lith  is 
considered  one  of  the  most  promising  geothermal  targets  that  encounter  many  hot  springs  with  a  good 
surface  temperature  upto  95  °C.  This  paper  aims  mainly  to  evaluate  the  geothermal  potential  of  the  main 
hot  spring  at  Wadi  Al-Lith  area  (Ain  A1  Harrah).  The  available  remote  sensing  images  are  analyzed  and  a 
number  of  2D  electric  resistivity  profiles  are  interpreted  to  delineate  the  surface  geological  lineaments 
and  the  subsurface  structural  elements  that  control  the  movement  of  the  thermal  water.  It  is  found  that 
the  main  surface  lineaments  are  structurally  oriented  along  NNE-SSW  and  NE-SW  directions  with  a 
frequency  percentage  of  52%  and  an  average  lineament  length  of  835  m.  Furthermore  the  subsurface 
structural  elements,  as  inferred  from  the  interpretation  of  the  geophysical  2D  electric  profiles,  have 
assumed  the  same  directions  beside  the  NW-SE  direction. 

The  characteristics  of  thermal  water  of  the  hot  spring  are  indicated  through  analyzing  the  major  and 
minor  elements  of  some  collected  water  samples.  Geo-thermometers  are  applied  to  estimate,  subsurface 
temperature,  heat  flow  and  discharge  enthalpy.  These  parameters  are  found  to  be  136  °C,  183  mW/M* 1 2 3  and 
219  kj/kg,  respectively.  An  estimate  of  the  geothermal  reserve  using  the  volumetric  method,  gave  total 
stored  heat  energy  of  1.713  x  1017  J  (rock  and  fluid)  and  a  geothermal  reserve  potential  of  26.99  MWt. 

It  appears  from  our  research  that  the  estimated  energy  is  quite  enough  for  operating  a  medium  scale 
power  plant  that  utilizes  low  boiling  point  fluid  (Kalina  Cycle)  for  limited  electricity  production,  beside 
other  low- temperature  applications  (district  heating,  green  houses  and  medical  therapy). 

©  2013  Elsevier  Ltd.  Open  access  under  CC  BY-NC-SA  license. 
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1.  Introduction 

Geothermal  power  is  among  the  renewable  energy  sources  which 
come  from  natural  resources  such  as  sunlight,  wind,  rain,  tides,  and 
biomass  which  are  renewable.  About  16%  of  global  final  energy 
consumption  comes  from  renewables,  with  10%  coming  from  tradi¬ 
tional  biomass,  which  is  mainly  used  for  heating,  and  3.4%  from 
hydroelectricity.  New  renewables  (small  hydro,  modern  biomass, 
wind,  solar,  geothermal,  and  biofuels)  accounted  for  another  3%  and 
are  growing  very  rapidly.  The  share  of  renewables  in  electricity 
generation  is  around  19%,  with  16%  of  global  electricity  coming  from 
hydroelectricity  and  3%  from  new  renewables  [1,2]. 

The  most  critical  factor  for  the  classification  of  geothermal  energy 
as  a  renewable  energy  scale  is  the  rate  of  energy  recharge.  In  the 
exploitation  of  natural  geothermal  systems,  energy  recharge  takes 
place  by  advection  of  thermal  water  on  the  same  time  scale  as 
production  from  the  resource.  This  justifies  our  classification  of 
geothermal  energy  as  a  renewable  energy  resource.  In  the  case  of 
hot,  dry  rocks,  and  some  of  the  hot  water  aquifers  in  sedimentary 
basins,  energy  recharge  is  only  by  thermal  conduction;  due  to  the  slow 
rate  of  the  latter  process,  however,  hot  dry  rocks  and  some  sedimen¬ 
tary  reservoirs  should  be  considered  as  finite  energy  resources  [3]. 

Geothermal  power  is  cost  effective,  reliable,  sustainable,  and 
environmentally  friendly,  but  has  historically  been  limited  to  areas 
near  tectonic  plate  boundaries.  Recent  technological  advances 
have  dramatically  expanded  the  range  and  size  of  viable  resources, 
especially  for  applications  such  as  home  heating,  opening  a 
potential  for  widespread  exploitation.  Geothermal  wells  release 
greenhouse  gases  trapped  deep  within  the  earth,  but  these 
emissions  are  much  lower  per  energy  unit  than  those  of  fossil 
fuels.  The  International  Geothermal  Association  (IGA)  has  reported 
that  10,715  MW  of  the  online  geothermal  power  in  24  countries, 
are  generating  67,246  GWh  of  electricity  in  2010.  This  represented 
a  20%  increase  in  geothermal  power  online  capacity  since  2005. 
IGA  expects  that  this  will  grow  to  18,500  MW  by  2015,  due  to  the 
large  number  of  projects  presently  under  consideration,  often  in 
areas  previously  assumed  to  have  little  exploitable  resource  [4,5]. 

In  Saudi  Arabia  there  is  a  big  gap  of  knowledge  and  information 
concerning  the  potential  of  geothermal  resources.  A  more  detailed 
and  much  more  comprehensive  work  is  needed  to  study  the 
geothermal  resources  in  Saudi  Arabia,  especially  those  encountered 
in  Southwestern  parts,  from  one  hand,  and  to  evaluate  their 
economic  reserves  for  possible  energy  production  and  domestic 
application,  from  the  other  hand.  Large  number  of  exploration 
methods  and  technologies  should  be  initiated  in  order  to  reach  these 
objectives.  Many  of  these  methods  are  in  current  use  and  have 
already  been  widely  experimented  in  other  sectors  of  research  [6]. 

The  geothermal  recourses  of  Saudi  Arabia  can  be  categorized  as 
the  following: 

1.  Low  enthalpy  resources  (sedimentary  aquifers):  represented  by 

deep-seated  aquifers  encountered  in  thick  sedimentary  basins 

in  the  eastern  part  of  the  Kingdom.  These  resources  are 


confined  and  geothermal  potential  is  represented  by  the 
normal  geothermal  gradient.  It  can  be  accessed  only  by  deep 
drilled  wells  (i.e.  oil  wells). 

2.  Medium  enthalpy  resources  (hot  springs):  encountered  along  the 
western  and  southwestern  coastal  parts  and  represented  by  the 
shallow  hot  springs  of  hot  surface  water.  These  resources  are 
unconfined  targets  with  direct  accesses  to  the  subsurface  hot 
anomalies  through  an  open  network  of  active  faults  and 
fractures  (structure  control). 

3.  High  enthalpy  resources  (basaltic  lavas ,  Harrats):  Saudi  Arabia 
has  approximately  80,000  km2  of  lava  fields,  known  as  Harrats. 
They  are  represented  by  volcanic  eruptions,  mainly  basaltic  in 
composition,  that  extend  along  the  coastal  part  of  the  Red  Sea 
at  the  western  of  Saudi  Arabia.  Harrats  of  Khaybar  and  Rahat 
are  believed  to  be  the  best  in  terms  of  high  heat  flow  and 
enthalpy  [6]. 

The  first  geothermal  work  done  at  Saudi  Arabia  was  that  of  El 
Dayel,  1988  [7].  He  focused  mainly  on  studying  the  hydro-chemical 
characteristics  of  the  hot  springs  at  Jizan  and  Al-Lith  areas.  Some 
water  samples  from  the  hot  springs  were  analyzed.  Other  inter¬ 
ested  work  regarding  the  petrography,  petrochemical  composition 
and  mineralogy  of  the  basaltic  volcanic  eruptions,  was  done  by 
Roobol  et  al.  and  Pint  et  al.  in  the  period  from  1992  to  2007.  A 
more  detailed  and  recent  work  is  done  by  Lashin  and  Al  Arifi,  2012 
[8].  They  focused  on  giving  a  preliminary  investigation  of  the 
geothermal  resources  encountered  in  the  southwestern  parts  of 
the  Saudi  Arabia  (Jizan  area),  based  on  the  interpretations  of  the 
Landsat  images  and  chemical  analyses. 

In  general,  a  few  studies  have  presented  the  characteristics  and 
the  potentiality  of  the  geothermal  resources  and  its  utility  for 
some  locations  along  the  western  coastal  area  of  Red  Sea  of  Saudi 
Arabia  [8-14].  This  work  presents  an  initial  evaluation  of  the 
geothermal  energy  potential  of  Wadi  Al-Lith  area,  and  provides  a 
basis  for  a  preliminary  assessment  of  geothermal  resources  in  the 
whole  region. 

2.  Geologic  setting 

2.1.  Rock  units 

Wadi  Al-Lith  is  represented  by  a  catchment  area  that  is 
considered  an  integral  part  of  the  Arabian  Shield.  It  extends 
from  the  western  coast  of  the  Red  Sea  to  high  Mountains  in  the 
east.  The  area  is  mainly  covered  by  metavolcanic  rocks,  metasedi¬ 
ments  and  late  Proterozoic  plutonic  rocks  (Fig.  1).  Four  major  rock 
units  are  found  covering  about  91%  of  Wadi  Al-Lith  catchment 
area.  These  rock  units  are  summarized,  from  youngest  to  oldest,  as 
follows: 

1-  Quaternary,  sand,  gravel  and  silt:  This  unit  represents  18%  of  the 
total  catchment  area  of  Wadi  Al-Lith.  Eolian  sand-dune  fields, 
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Fig.  1.  Geologic  map  of  Wadi  Al-Lith  area  showing  Ain  Al  Harrah  hot  spring. 


mobile  sheet  sand  and  silt  deposits  cover  a  large  proportion  of 
the  coastal  plain  in  the  southern  part  of  the  study  area.  Barrier, 
longitudinal,  and  barchans  dunes  are  present  locally,  and  most 
of  the  area  is  covered  by  sand  deposits. 

2-  Late-  to  post-tectonic  granitic  rocks :  These  rocks  are  part  of  the 
Proterozoic  intrusive  rocks  that  cover  most  of  the  western  part 
of  the  study  area.  They  represent  13%  of  the  total  basin  area  of 
Wadi  Al-Lith.  The  plutonic  rocks  range  in  composition  from 


serpentinite  to  syenite,  although  diorite,  tonalite,  granodiorite, 
and  monzogranite  predominate. 

3-  Lith  suite,  Khasrah  complex,  diorite  and  gabbro :  These  rocks 
contain  mafic  to  intermediate  plutonic  rocks  that  crops  out  in  a 
wide  northeast  trending  belt  in  the  western  part  of  Wadi  Al- 
Lith  basin.  It  covers  almost  37%  from  Wadi  Al-Lith  basin  area. 

4-  Baish  and  Baha  groups :  Many  rocks  assigned  to  the  undivided 
Baish  and  Bahah  groups  crops  out  in  the  eastern  part  of  Wadi 
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Al-Lith  basin  as  “basalt  and  dacite”,  and  “biotite-hornblende 
schist,  and  amphibolite”  covering  almost  23%  from  Wadi  Al-Lith 
basin  area. 

Other  rock  units  are  found  constituting  about  9%  of  wadi  Al- 
Lith  basin.  The  most  important  of  these  rocks  are  the  metagabbro, 
gabbro  dikes,  syentectonic  granitic  rock,  granodiorite,  granite  and 
amphibolite,  as  shown  on  the  geological  map  (Fig.  1 ). 


The  main  structural  elements  which  are  prevailing  along  the 
western  coastal  part  (Red  Sea)  of  the  Saudi  Arabia,  including  the 
study  area,  can  be  classified  as: 

•  Precambrian  northwest-southeast  to  east-west  compression 
that  are  reactivated  by  tensional  stress  field  during  the  Tertiary 
evolution  of  the  Red  Sea  rift. 

•  Tertiary  northeast-southwest  tension. 
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Fig.  2.  Digital  elevation  model  of  Wadi  Al-Lith  basin. 
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These  structural  elements  are  believed  to  be  the  main  channels 
and  pathways  through  which  the  thermal  water  rises  to  the  surface. 

2.2.  Topography 

Fig.  2  exhibits  the  digital  elevation  model  of  Wadi  Al-Lith  basin, 
which  shows  a  wide  variation  in  the  topography  of  the  area. 
A  maximum  depth  of  2600  m  (above  sea  level)  is  recorded  at 
the  highly  elevated  mountainous  area  in  the  northern  part  of  Wadi 


Al-Lith  basin,  while  an  elevation  of  less  than  1  m  is  recorded  for 
the  sea  warded  coastal  plains  in  the  southern  part  of  the  area. 

According  to  its  surface  elevation  from  the  sea  level,  Wadi 
Al-Lith  basin  can  be  divided  into  three  zones,  as  follows: 

1-  A  high  mountainous  upstream  zone:  the  elevation  ranges 
from  500  m  to  2600  m  above  the  sea  level  and  is  characterized 
by  high  degrees  of  slopes  and  rugged  area.  It  includes  a  number 
of  tributaries  and  sub-wadies  such  as;  Wadi  Birayn,  Wadi 
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Fig.  3.  Landsat  RGB  image  of  Wadi  Al-Lith  catchment  area. 
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Salibah  and  Wadi  Maqsa  which  join  the  mainstream  of  Wadi 
Al-Lith. 

2-  An  intermediate  zone  of  small  hills  and  mountains:  the  eleva¬ 
tion  of  this  zone  ranges  from  100  m  to  500  m  above  sea  level.  It 
represents  the  midstream  of  the  Wadi  Al-Lith  and  includes  the 
main  tributaries  and  wadies  sloping  from  up-stream.  It  is 
characterized  by  moderate  gradient  and  encounters  the  major 
waterways  of  the  Wadi  Al-Lith. 

3-  The  coastal  flood  plain  down-stream  zone:  this  zone  contains 
alluvial  deposits  with  an  elevation  range  from  100  m  to  nearly 
zero  at  the  shoreline  of  the  Red  Sea. 


3.  Methodology 

The  utilized  methods  in  this  work  are  mainly  applied  to  detect 
the  following  issues: 

1-  The  surface  and  subsurface  structural  elements  that  control  the 
movement  of  the  ascending  thermal  water. 

2-  The  petro-thermal  characteristics  of  thermal  water. 

3-  The  possible  geothermal  energy  potential  and  its  reserves. 


The  analysis  of  surface  and  subsurface  structural  elements  are 
conducted  through  analyzing  the  available  remote  sensing  images 
and  conducting  a  number  of  geophysical  2D  electric  profiles. 


3.1.  Landsat  RGB  images 

The  remote  sensing  imagery  is  useful  in  mapping  lithology  and 
offers  an  efficient  tool  for  detecting  geological  structures,  such  as 
lineaments.  In  the  present  study,  extraction  of  lineaments  is  enhanced 
through  analyzing  the  available  Landsat  RGB  images  with  30  m 
resolutions  (Fig.  3).  PCI-Geomatica  software  is  used  as  a  powerful  tool 
for  extracting  the  surface  geological  structures  in  the  study  area. 

3.2.2.  Gridding  and  mapping  oflineation 

The  lineation  data  that  are  extracted  from  the  analysis  of 
Landsat  images  are  feed  to  specific  software  (RockWorks  Utilities) 
to  create  a  grid  file  that  models  all  aspects  of  lineation  frequency, 
cumulative  lengths,  and/or  intersections,  with  a  variety  of  weight¬ 
ing  options.  Three  main  maps  were  constructed: 

1.  Lineament  intersection  map :  It  illustrates  the  intersection  grid 
map  in  which  each  cell  represents  the  number  of  intersections 
of  the  lineations  that  occur  within  the  cell. 

2.  Lineament  length  map :  It  is  a  length  grid  map,  in  which  each  cell 
represents  the  total  added  lengths  of  all  of  the  lineations  that 
pass  through  this  cell  (total  length  of  each  lineation  is 
included). 

3.  Lineament  frequency  map:  It  is  a  lineation  frequency  grid  map 
that  includes  the  number  of  lineations  that  lie  in  the  cell. 

On  the  other  hand,  a  2D  geophysical  electric  resistivity  method, 
is  applied  to  investigate  the  subsurface  orientation  of  the 


Table  1 

The  chemical  analyses  of  Ain  Al  Harrah  hot  spring. 


Locality 

Wadi  Al-Lith 

Sampling  date 

28  May  2011 

TDS 

2180  ppm 

Name  &  No. 

Ain  Al  Harrah 

Type  of  aquifer 

Hot  spring 

PH 

7.8 

Longitude 

40°28/16.4// 

Temperature 

<90  °C 

EC 

3633  pS  cm  1 

Latitude 

20°27/41// 

Elevation 

167.0  m 

Chemical  composition  (mg/1) 

Si02 

Na+  I<+ 

Ca++ 

Mg++ 

i 

i 

o 

C/5 

nh4 

HC03“ 

cr 

F" 

72.69 

403.5  25.2 

260 

0.461 

225.4 

* 

18.481 

845 

2.53 

Chemical  composition  (10  3  mg/1) 

Li 

Al  Zn 

Fe 

Cu 

Mn 

Ni 

Co 

Ba 

B 

90.5 

6.9  3.39 

647 

3.66 

1.35 

3.8 

0.393 

80.9 

378 

Chemical  composition  (ppb) 

Rb 

Cr  V 

Mo 

Sr 

U 

Th 

Not  detected 

68.9 

87.4  28.9 

34.3 

4160 

0.216 

0.0217 

Over  range 

Table  2 

The  chemical  analyses  of  water  well  close  to  the  hot  spring. 

Locality 

Wadi  Al-Lith 

Sampling  date 

28  May  2011 

TDS 

1127  ppm 

Name  &  No. 

28 

Type  of  aquifer 

Water  well 

PH 

7.5 

Longitude 

20°17/32// 

Temperature 

34  °C 

EC 

1879  pS  cm-1 

Latitude 

40°28/36.1// 

Elevation 

60.0  m 

Chemical  composition  (mg/1) 

Si02 

Na+ 

K+ 

Ca++ 

Mg++ 

i 

i 

O 

C/5 

NH. 

4 

HC03“ 

cr 

F“ 

39.28 

212.6 

9.34 

103 

37.8 

340 

38.872 

185.9 

0.61 

Chemical  composition  (10  3  mg/1) 

Li 

Al 

Zn 

Fe 

Cu 

Mn 

Ni 

Co 

Ba 

B 

1.68 

9.21 

2.93 

270 

1.91 

0.638 

4.06 

0.294 

27.4 

505 

Chemical  composition  (ppb) 

Rb 

Cr 

V 

Mo 

Sr 

U 

Th 

Not  detected 

0.587 

56.8 

41 

23.2 

1630 

2.6 

Over  range 
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Fig.  4.  Interpretation  of  the  lineaments  of  Wadi  Al-Lith  area  as  extracted  from  the  analysis  of  the  Landsat  GRB  image. 
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Fig.  5.  Rose  diagram  showing  the  structure  lineaments  of  (A)  Wadi  Al-Lith  area  and  (B)  Ain  Al  Harrah  hot  spring  area. 
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Table  3 

The  directional  trends  of  the  extracted  lineaments  in  the  Wadi  Al-Lith  area. 


Lineaments  structure  trend 

Degrees  range 

Frequencies 

Percent 

Avg.  length  (m) 

SD 

N-S 

(0-11.25),  (168.75-191.25)  and  (348.75-360) 

1644 

16.2 

759 

569 

NNE-SSW 

(11.25-33.75)  and  (191.25-213.75) 

2297 

22.6 

745 

488 

NE-SW 

(33.75-56.25)  and  (213.75-236.25) 

2230 

21.9 

764 

535 

ENE-WSW 

(56.25-78.75)  and  (236.25-258.75) 

1683 

16.5 

757 

619 

E-W 

(78.75-101.25)  and  (258.75-281.25) 

839 

8.2 

777 

699 

ESE-WNW 

(101.25-123.75)  and  (281.25-303.75) 

329 

3.2 

1287 

1153 

SE-NW 

(123.75-146.25)  and  (303.75-326.25) 

339 

3.3 

1095 

827 

SSE-NNW 

(146.25-168.75)  and  (326.25-348.75) 

815 

8.0 

968 

676 

Table  4 

The  directional  trends  of  the  extracted  lineaments  of  Ain  Al  Harrah  area. 

Lineaments  structure  trend 

Degrees  range 

Frequencies 

Percent 

Avg.  length  (m) 

SD 

N-S 

(0-11.25),  (168.75-191.25)  and  (348.75-360) 

99 

13.3 

803 

528 

NNE-SSW 

(11.25-33.75)  and  (191.25-213.75) 

193 

25.9 

822 

613 

NE-SW 

(33.75-56.25)  and  (213.75-236.25) 

195 

26.1 

847 

644 

ENE-WSW 

(56.25-78.75)  and  (236.25-258.75) 

97 

13.0 

687 

435 

E-W 

(78.75-101.25)  and  (258.75-281.25) 

33 

4.4 

922 

757 

ESE-WNW 

(101.25-123.75)  and  (281.25-303.75) 

24 

3.2 

164 

83 

SE-NW 

(123.75-146.25)  and  (303.75-326.25) 

61 

8.2 

154 

54 

SSE-NNW 

(146.25-168.75)  and  (326.25-348.75) 

44 

5.9 

214 

89 
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Fig.  6.  Lineation  intersection  grid  map  of  Wadi  Al-Lith  basin  including  Ain  Al  Harrah  area. 
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Fig.  7.  Lineation  length  grid  map  of  Wadi  Al-Lith  basin  including  Ain  Al  Harrah  area. 


geothermal  targets  and  detecting  their  possible  lateral  extension 
as  well  as  the  structural  elements  which  may  control  the  subsur¬ 
face  movement  of  the  arising  geothermal  water. 

3.2.  Geophysical  surveys 

The  geothermal  potential  can  be  explored  and  investigated 
using  a  variety  of  geophysical  and  geochemical  (geo-thermo- 
meters)  techniques  which  are  normally  used  in  geo-technical 
and  geological  investigation  as  well  as  in  the  exploration  of 
hydrocarbons  in  the  oil  industry  [14-22].  Such  techniques  provide 
the  necessary  information  for  constructing  a  realistic  model  of  the 
geothermal  system  and  assessing  the  potential  of  the  resource.  The 
electric  resistivity  is  related  to  various  geological  parameters  such 
as  the  mineral  and  fluid  content,  porosity  and  degree  of  water 
saturation  in  the  rock  [23]. 

A  number  of  2D  electric  profiles  are  performed  at  the  location 
of  the  hot  springs  in  Wadi  Al-Lith  area.  The  Syscal-Rl  system  (IRIS 
instruments)  with  72  multi-electrodes  system,  with  4  x  18  elec¬ 
trode  cables  is  used  for  conducting  the  field  surveys  using  different 
electrode  spacing  (1  m,  3  m,  5  m  and  10  m)  arranged  along  a 


profile.  Wenner-Schlumberger  arrangement  is  used  and  the  sur¬ 
veyed  sections  are  taken  very  close  to  and  crossing  the  hot  spring. 
RES2DINV  software  is  used  for  processing,  correcting  and  inter¬ 
pretation  of  electric  resistivity  data.  It  can  operate  large  datasets 
collected  by  a  large  number  of  electrodes  and  can  account  for  the 
topographical  effect  along  the  survey  lines.  Finally,  it  generates 
two-dimensional  resistivity  model  for  subsurface  from  electrical 
imaging  surveys  [23]. 


3.3.  Petro-thermal  characteristics  of  thermal  water 

The  petro-thermal  characteristics  of  thermal  water  of  Ain  Al 
Harrah  hot  spring  at  Wadi  Al-Lith  area  are  mainly  enhanced  through 
analyzing  water  samples  collected  from  the  hot  spring  and  a  very 
close  water  well.  These  water  samples  are  analyzed  using  the 
ICP-MS  spectrometer.  The  analysis  was  performed  by  Perkin-Elmer 
Sciex  instruments  multi-element  ICP-MS  spectrometer  (type 
ELAN9000)  equipped  with  standard  torch,  cross  flow  nebulizer 
and  Ni  sampler  and  skimmer  cones.  Tables  1  and  2  show  the 
chemical  analyses  (major,  minor  and  trace  elements)  of  two  water 
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Fig.  8.  Lineation  frequency  grid  map  of  Wadi  Al-Lith  basin  including  Ain  Al  Harrah  area. 


samples  taken  from  Ain  Al  Harrah  hot  spring  and  the  neighboring 
water  well. 


3.3 A.  Geo-thermometers 

Geo-thermometer  is  a  chemical  substance  or  isotope  in  water 
that  may  be  used  to  predict  the  subsurface  temperature  in 
geothermal  systems.  The  basic  assumption  behind  using  the  geo¬ 
thermometers  is  the  chemical  equilibrium  between  the  substance 
(isotope)  and  mineral  (s)  in  the  reservoir.  Theoretically,  any  cation 
ratio  and  any  aqueous  species  concentration  can  be  used  as  a  geo¬ 
thermometer  as  long  as  equilibrium  prevails.  A  temperature 
equation  for  a  geo-thermometer  is  a  temperature  equation  for  a 
specific  equilibrium  constant  referring  to  a  specific  mineral- 
solution  reaction  [24-29]. 

A  selected  number  of  geo-thermometers  (Quartz,  Chalcedony, 
Na,  K,  Ca  geo-thermometers  and  combinations  of  them;  i.e.  Na-Ca 
and  Na-K-Ca)  are  applied  for  determining  the  subsurface  forma¬ 
tion  temperature,  discharge  enthalpy  and  heat  flow.  These  geo¬ 
thermometers  have  a  low  to  medium  temperature  range  of  25- 
250  °C  that  fit  with  the  thermal  properties  of  the  hot  spring  under 


study.  Eqs.  (l)-(7)  demonstrate  the  mathematical  expressions  of 


these  geo-thermometers  [24-29]: 

(a)  Silica  geo-thermometers 
Quartz 

Ti(°C)  =  5.i9-”gQz-273-15  (Range  of  °C  :  25-250)  (1) 

r2(0C)  =  42.198  +  0.2883  lQz-3. 6686  x  10 ~4Qz2  +  3.1665 

x10“7Qz3  +  77.034LogQz  (Range  of  °C  :  25-250)  (2) 

Chalcedony 

Ti(°Q=  4  69!°32gQz-273.15  (Range  of  °C  :  25-250)  (3) 

1112 

T2CQ  =  4  91-LogQz~273  15  <Range  °f  °C  1  25-250)  (4} 


(b)  Na-K-Ca  geo-thermometers 
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Na-I< 

^i(°C)=  o.993+Log(Na/K) -273.15  (Range  of  °C  :  25-250)  (5) 

Na-C 

T'(°C>  =  3.08-LogtNa/ca05)  —273.1 5  (Range  of  X  :  25-250)  (6) 

Na-K-Ca 

Ti(°C)  = - 1647  - 273.15  (Range  of  °C  :  0-350) 

Log(Na//<')+/?Log(Ca°'5/Na)+2.24  \  B  ) 

(7) 

3.32.  Ternary  diagrams 

Three  ternary  (C1-S04-HC03  and  Na+K-Ca-Mg)  and  Giggen- 
bach  (Na-K-Mg)  diagrams  are  constructed  for  Ain  Al  Harrah  hot 
spring  to  classify  its  thermal  water  (based  on  major  anions  and 
cations)  and  to  indicate  the  prevailing  subsurface  thermal  condi¬ 
tions  [30,31]. 


3.4.  Geothermal  reserve  estimation 

Geothermal  reserves  are  defined  as  quantities  of  thermal 
energy  which  are  expected  to  be  recovered  from  known  reservoirs 
from  a  given  date  forward.  The  reserve  is  that  part  of  the  resources 
which  could  be  extracted  economically  and  legally  at  present  and 
that  are  known  and  characterized  by  drilling  or  by  geochemical, 
geophysical  and  geological  evidence  [32,33]. 

The  most  common  approach  for  geothermal  reserve  estimation 
is  the  volumetric  method.  This  method  is  mainly  utilized  in  the 
later  stages  of  geothermal  exploration  and  after  that  in  the  early 
stages  of  field  development  to  justify  drilling  and  commitment  for 
a  specified  power  plant  size,  although  may  be  used  if  good  surface 
indications  and  geophysical  surveys  are  available  [34].  It  can  be 
better  applied  during  this  stage  than  numerical  modeling  which 
requires  a  significant  number  of  wells  and  production  history  to  be 
considered  reliable  [35]. 

In  general,  the  reservoir  parameters  necessary  for  geothermal- 
reserve  estimation  are  gathered  from  the  analyses  of  the  tem¬ 
perature  data  and  from  other  geophysical  and  geochemical 
interpretations.  A  computer  program  (RESPAR,  ICEBOX  software) 
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Fig.  9.  Accumulation  grid  map  of  the  intersection,  length  and  frequency  of  the  extratced  lineaments  of  Wadi  Al-Lith  basin. 
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is  used  to  estimate  the  different  reservoir  properties  providing 
that  the  temperature,  porosity  and  the  reservoir  thickness 
parameters  are  known.  The  total  stored  energy  of  the  geothermal 
reservoir  can  be  regarded  as  the  sum  of  energy  stored  in  the  rock 
(Er)  and  the  energy  represented  by  the  fluid  stored  in  the  pore 
spaces  (Ef).  It  can  be  estimated  using  the  following  equation: 

Et  =  Er  +  Ef  =  VO-<f>)prCr(Ti-T0)  +  V<p  pwCw(T[-T0)  (8) 

Where  Et  is  the  total  thermal  energy  (J)  in  the  rock  (Er)  and  fluid 
(Ef);  cf)  is  the  reservoir  porosity  (%);  V  is  the  reservoir  volume  (m3); 
/9r>w  are  the  densities  of  rock  and  water  (kg/m3);  Cr>w  are  the  heat 
capacities  of  rock  and  water  (J/kg  °C)  and  TQ  are  the  initial 
reservoir  and  the  reference  temperatures  (°C). 

This  stored  thermal  energy  can  be  converted  into  power 
potential  using  the  following  equation: 

Et-  x  RF  x  CF 

Power  Potential  (MWt)  =  — — — — —  (9) 

PL  x  Lr 

where  RF  is  the  recovery  factor,  CE  is  the  conversion  efficiency,  PL 
is  the  geothermal  plant  life  in  years  and  LF  is  the  load  factor. 


Fig.  10.  Base  map  of  Ain  Al  Harrah  hot  spring,  showing  the  location  of  the  2D 
electric  profiles. 


4.  Results  and  discussion 

4.2.  Lineaments  extraction  and  analysis 

The  analysis  of  lineaments  has  clarified  that  about  10,000 
elements  are  detected  in  Wadi  Al-Lith  basin  (Fig.  4).  RockWorks 
software  is  mainly  utilized  for  drawing  the  rose  diagram  for 
detecting  the  general  trends  of  these  lineaments  in  Wadi  Al-Lith 
basin  and  more  specifically  in  Ain  Al  Harrah  area  (Fig.  5).  The 
diagram  illustrates  that  the  prevailing  structural  elements  are 
nearly  the  same  and  almost  all  lineaments  are  oriented  and 
trended  from  0°  to  90°  (or  180°-270°)  with  mean  trend  vector  of 
33°  (or  213°).  Few  lineaments  are  oriented  along  a  northwest- 
southeast  direction  in  Ain  Al  Harrah  hot  spring  area  (Fig.  5b). 
Tables  3  and  4  summarize  the  main  structural  trends  of  these 
lineaments,  as  well  as  their  percent,  lengths  and  frequency. 

Figs.  6-8  show  the  constructed  intersection,  length  and  fre¬ 
quency  property  maps  for  all  the  extracted  lineaments  in  Wadi  Al- 
Lith  area.  An  accumulation  grid  map  for  these  three  properties  is 
drawn,  to  illustrate  their  combined  effect  (Fig.  9).  This  map  shows 
how  much  these  three  components  behave  and  affect  the  study 
area.  The  weighting  factor  or  influence  of  each  component  in  the 
accumulation  grid  map  was  chosen  to  be,  lineation  intersection 
(3),  lineation  length  (2)  and  lineation  frequency  (1). 

Interpretation  of  these  figures  shows  that,  no  more  than  three 
lineament  elements  are  acting  together  in  the  study  area.  This 
assigns  middle  to  low  influence  of  these  elements  in  Wadi  Al-Lith 
catchment  area.  Outside  Wadi  Al-Lith  area,  on  the  other  hand,  a 
high  combined  influence  of  these  structural  lineaments  is 
recorded  in  the  order  of  four  elements  (Fig.  9).  The  result  that 
can  be  concluded  from  this  lineament  analysis  is  that  the  surface 
structural  elements  are  not  the  factor  that  control  the  movement 
of  the  deep  arising  thermal  water  of  the  main  hot  spring  in  Wadi 
Al-Lith  area,  and  that  there  might  be  another  subsurface  elements 
that  could  carry  out  this  geothermal  water  to  the  surface. 

4.2.  Interpretation  of  geophysical  profiles 

Out  of  the  nine  conducted  2D  electric  resistivity  sections,  three 
are  selected  to  demonstrate  the  prevailing  structural  elements. 
These  profiles  are  specifically  conducted  in  the  area  of  the  hot 
spring  (Fig.  10).  The  results  obtained  from  these  profiles  are  mainly 
concerned  with  the  pathways  of  geothermal  water,  detecting  feed 
zones,  the  influencing  structural  elements  and  the  possible  facture 
system. 

Fig.  11  exhibits  the  interpreted  2D  electric  section  (Line  I)  that  is 
extended  along  an  east-west  direction.  It  shows  a  tilted  thermal 
feed  zone  at  the  middle  of  the  section  oriented  along  a  very 
obvious  fault  zone.  The  strike  of  this  fault  is  found  to  be  along  a 
northwest-southeast  direction  which  is  considered  the  main 
prevailing  structure  in  the  area  (Red  Sea  Rift).  This  structure  trend 
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Fig.  11.  2D  interpreted  resistivity  profile  (line  I)  extending  East-West,  Wadi  Al-Lith  area. 
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Depth  Iteration  7  fibs,  error  =  2.9  7. 
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Fig.  12.  2D  interpreted  resistivity  profile  (line  H)  extending  East-West,  Wadi  Al-Lith  area. 
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Fig.  13.  2D  interpreted  resistivity  profile  (line  B)  extending  North-South,  Wadi  Al-Lith  area. 


Table  5 

Petro-thermal  parameters  of  Ain  Al  Harrah  hot  spring,  Wadi  Al-Lith  area. 


Hot  spring 

Subsurface  temp.  (°C) 

Discharge  enthalpy  (kj/kg) 

Heat  flow  (mW/M2) 

Ain  Al  Harrah 

136 

219 

183 

0  Cl  100 

B 

Na+K 

Fig.  14.  Type  of  thermal  fluids  of  Ain  Al  Harrah  hot  spring  as  inferred  from;  (A)  C1-S04-HC03  and  (B)  NA+K-Ca-Mg  ternary  diagrams. 


is  opposite  to  another  oriented  northeast-southwest  geothermal 
feed  zone,  detected  at  the  western  end  of  the  profile  H  between 
electrodes  80-90  (Fig.  12). 


Fig.  13,  on  the  other  hand,  represents  another  2D  electric  profile 
extending  from  north  to  south  (Line  B).  A  suggested  fault  system  is 
recognized  at  the  middle  of  the  section.  This  may  be  interpreted  as  a 
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result  of  the  well  known  northeast-southwest  trending  faults  that 
influence  the  western  coastal  parts  of  the  Red  Sea  area.  In  general, 
the  feed  zones  of  thermal  waters  are  indicated  by  lower  resistivity 
values  and  the  structural  control  of  the  thermal  water  is  obvious. 

43.  Characteristics  of  thermal  water 

The  petro-thermal  parameters  which  are  indicated  from  the 
geo-thermometers  are  listed  in  Table  5.  It  appears  clearly  that  the 
Ain  Al  Harrah  hot  spring  attains  good  and  high  geothermal 
parameters.  The  weighted  values  of  the  estimated  subsurface 
formation  temperature,  discharge  enthalpy  and  heat  flow  are 
found  to  be  136  °C,  219  kj/kg  and  183  mW/M2,  respectively. 

Figs.  14  and  15  represent  the  constructed  ternary  (C1-S04-HC03 
and  Na+K-Ca-Mg)  and  Giggenbach  (Na-K-Mg)  diagrams  for  Ain  Al 
Harrah  hot  spring  and  one  neighboring  water  well.  These  dia¬ 
grams  are  used  to  indicate  the  type  of  the  geothermal  water  and 
the  subsurface  geothermal  condition  at  which  the  dissolved  ions 
of  the  surface  ascending  thermal  fluids  are  originated. 

The  ternary  diagrams  illustrate  that  the  water  type  of  the  hot 
spring  is  completely  different  than  that  of  the  neighboring  water 
wells.  The  increase  of  the  concentration  of  the  chlorine  anion  (mg/1) 
on  the  expanse  of  the  sulfate  and  carbonate  inions  is  very  clear  on 
the  hot  spring  (Fig.  14a).  On  the  other  hand,  a  considerable  increase 
of  the  Mg  cation  is  indicated  for  the  water  well  sample  on  the 
expanse  of  the  Na+K  and  Ca  cations  (Fig.  14b).  A  mature  water  type, 
most  probably  mixed  with  surface  water  is  indicated  to  the  hot 
spring,  while  a  volcanic  unmixed  water  type  is  given  for  the  water 
well.  Giggenbach  (Na-K-Mg)  diagram  shows  that  Ain  Al  Harrah  hot 
spring  assigns  higher  geothermal  parameters  than  the  water  well. 
The  cluster  of  this  hot  spring  is  located  along  the  150  °C  Mg-I<  and 
the  175  °C-225  °C  Na-K  iso-thermal  lines,  while  that  of  the  water 
well  is  located  very  close  to  the  Mg  corner  along  iso-thermal  lines  of 
less  than  50  °C  assuming  very  low  geothermal  regime  (Fig.  15). 

Another  unique  difference  in  water  type,  is  the  increase  of  the 
concentration  of  some  radioactive  elements  in  the  hot  spring  like 
rubidium  (Rb,  68.9  ppm)  and  strontium  (Sr,  4160  ppm)  more  than 
those  in  the  water  well  (Rb,  0.587  ppm  and  Sr,  4160  ppm). 

4.4.  Geothermal  energy  potential 

In  this  study,  geothermal-reserve  estimation  is  carried  out  for 
evaluating  the  possible  geothermal  capacities  of  Ain  Al  Harrah  hot 


Table  6 

Reservoir  parameters  used  for  geothermal  reserve  estimation. 


Property 

Unit 

Value 

Reservoir 

Volume 

km3 

1.28 

Porosity 

% 

4 

Density 

Rock 

kg/m3 

2800 

Water 

kg/m3 

980 

Heat  capacity 

Rock 

J/kg  °C 

1000 

Water 

J/kg  °C 

4200 

Temperature 

Initial 

°c 

95 

Reference 

°c 

45 

spring.  Reserve  estimation  is  based  on  the  above  mentioned  Eqs. 

(8)  and  (9),  assuming  the  following  parameters: 

•  A  reservoir  area  of  4.8  km2  (1.2  x  4.0  km2)  and  reservoir  volume 
of  1.28  km3  providing  a  reservoir  thickness  of  250  m  as  con¬ 
cluded  from  the  interpretation  of  resistivity  data. 

•  An  initial  temperature  of  95  °C  for  the  geothermal  reservoir. 

•  A  reference  temperature  of  45  °C  is  suggested  (due  to  the 
higher  average  annual  surface  temperature  in  Saudi  Arabia). 

•  Pore  volume  range  of  3-5  p.u.  for  the  fractured  granitic  rocks 
(average  value  of  4  p.u.). 

•  Density  range  of  2.7-2.90  g/cc  (average  value  of  2.80  g/cc). 

The  different  reservoir  parameters  that  are  derived  from  the 
analyses  of  temperature  and  geo-thermometer  data  and  can  be  used 
in  thermal  energy  estimation  are  represented  in  Table  6.  According  to 
these  parameters  and  by  applying  in  Eq.  (8),  the  total  stored  heat 
energy  of  Ain  Al  Harrah  geothermal  area  is  found  to  be  1.713  x  1017  J 
(1.34  x  1016J  for  thermal  fluids  and  1.57  x  1017  J  for  reservoir  rock). 

The  total  stored  heat  energy  of  Ain  Al  Harrah  hot  spring  can  be 
converted  into  power  energy  units  (megawatt  thermal  using  Eq. 

(9) )  for  possible  energy  production  and  other  low  geothermal 
utilizations.  Assuming  recovery  factor  of  0.2  (normal  factor  for 
geothermal  reservoirs),  turbine  conversion  efficiency  of  0.47, 
(initial  temperature/reservoir  temperature)  plant  life  of  20  years 
and  load  factor  of  0.95,  the  thermal  power  potential  of  Ain  Al 
Harrah  geothermal  reservoir  is  found  to  be  26.99  MWt.  This 
geothermal  reserve  is  considered  a  very  promised  and  can  be 
used  for  providing  the  Al  Lith  area  with  electricity  (clean  energy) 
on  a  long  term  basis.  We  strongly  recommend  exploiting  this 
geothermal  system  to  generate  electricity  for  the  first  time  in 
Saudi  Arabia  using  a  specific  power  plant  that  operates  thermal 
fluid  that  boils  at  low  temperature  (e.g.  Kalina  Cycle). 

4.4.1.  Utilization  of  geothermal  energy  in  Saudi  Arabia  and  economic 
and  social  impact 

The  estimated  geothermal  energy  (  MWt)  is  quite  enough  for 
possible  electricity  production  in  Saudi  Arabia.  Kalina  power  plants 
can  be  used  due  to  the  lower  boiling  point  (70  °C)  of  their  fluids. 
Having  a  higher  temperature  range  than  the  boiling  point  of  the 
Kalina  cycle,  the  thermal  water  of  Ain  Al  Harrah  hot  springs  can  be 
economically  used  for  power  generation.  However,  this  thermal 
water  can  be  utilized  directly  in  other  low  forms  of  applications 
(district  heating,  fish  farming,  agricultural  applications,  greenhouses, 
etc.).  Some  of  these  applications  are  already  now  in  use  (Gouba  area, 
central  Saudi  Arabia),  where  a  number  of  geothermal-based  natural 
therapy,  medical  and  swimming  pools  are  constructed. 

The  social  factor  of  using  the  geothermal  energy  implies  from 
the  fact  that  most  of  the  geothermal  resources  are  mainly  located 
in  the  extreme  western  and  south-western  borders  of  Saudi  Arabia 
along  the  coast  of  the  Red  Sea.  These  areas  besides  being  far  are 
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topographically  complicated,  not  so  easy  to  reach  and  need  more 
developments  in  its  infrastructures  and  the  level  of  the  residence 
Bedouins  people.  The  geothermal  based-energy  can  help  in  con¬ 
structing  new  modern  communities  fed  with  the  geothermal- 
based  energy,  where  no  high  voltage  transmission  lines  are 
available  in  the  vicinity  and  where  it  would  be  too  expensive  to 
connect  to  the  national  electric  grid. 

5.  Conclusions 

The  following  conclusions  may  be  drawn  according  to  the 
present  study: 

1.  Wadi  Al-Lith  area  contains  one  of  the  most  promised  geother¬ 
mal  systems  in  the  western  of  Saudi  Arabia.  It  is  located  near 
the  coast  of  the  Red  Sea  and  is  strongly  affected  by  its  tectonic 
activity  and  structural  regimes. 

2.  Ain  Al  Harrah  hot  spring  is  considered  the  main  geothermal 
target  in  Wadi  Al-Lith  area  that  attains  high  surface  tempera¬ 
ture  upto  95  °C  and  good  petro-thermal  characteristics  of 
136  °C  (reservoir  temperature),  219  kj/kg  (discharge  enthalpy) 
and  183  mW/M2  (heat  flow)  as  indicated  from  the  geo¬ 
thermometric  analysis. 

3.  The  analysis  of  Landsat  RGB  images  (lineaments  extraction)  and 
the  interpretation  of  the  electric  geophysical  profiles  have 
clarified  the  presence  of  many  surface  and  subsurface  struc¬ 
tural  elements.  Feed  zones  of  thermal  waters  are  indicated  by 
lower  resistivity  values  and  the  structural  control  of  the 
thermal  water  is  obvious.  Majority  of  the  encountered 
geothermal  feed  zones  are  tilted  and  aligned  along  a  north- 
east-southwest  (transform)  direction  and  little  bit  to  a  north¬ 
west-southeast  one  (Red  Sea  Rift). 

4.  Almost  all  surface  structural  lineaments  that  control  the  move¬ 
ment  of  the  geothermal  water  (Ain  Al  Harrah  hot  spring)  are 
found  following  the  main  structural  trends  of  Wadi  Al- 
Lith  area. 

5.  The  analysis  of  ternary  and  Giggenbach  diagrams  shows  that 
the  water  type  of  geothermal  system  is  completely  different 
than  the  neighboring  domestic  water  aquifers.  A  mature  mixed 
water  type  is  assigned  for  Ain  Al  Harrah  hot,  while  a  volcanic 
unmixed  water  model  is  given  for  the  neighboring  water  wells. 

6.  The  stored  heat  energy  is  found  to  be  1.34xl016J  for  the 
thermal  fluids  and  1.57  x  1017J  for  the  hot  rock. 

7.  A  geothermal  power  potential  of  26.99  MWt  is  estimated  for 
Ain  Al  Harrah  hot  spring  which  is  considered  good  and 
adequate  for  economic  investment. 

8.  A  medium  scale  power  plant  that  utilizes  fluids  of  low  boiling 
point  is  recommended  for  possible  electricity  generation  for 
the  first  time  in  Saudi  Arabia. 

9.  A  very  strong  and  positive  economic-social  impact  will  be 
obtained  upon  applying  such  geothermal  technology  in  the 
study  area  or  in  other  remote  areas  of  similar  geothermal 
systems.  It  will  enhance  the  infrastructure  and  build  up  a 
new  geothermal-feed  communities  based  on  clean  renewable 
energy. 
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